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abstract — We use the scattering matrix approach to derive generalized Bardeen-like formulae 
for the conductances between the contacts of a phase-coherent multiprobe conductor and a tunneling 
tip which probes its surface. These conductances are proportional to local partial densities of states, 
called injectivities and emissivities. The current and the current fluctuations measured at the tip 
are related to an effective local non- equilibrium distribution function. This distribution function 
contains the quantum-mechanical phase-coherence of the charge carriers in the conductor and is 
given as products of injectivities and the Fermi distribution functions in the electron reservoirs. 
The results are illustrated for measurements on ballistic conductors with barriers and for diffusive 
conductors. 



Introduction — Over the last 15 years the scanning tunneling microscope (STM) has developed into a standard 
tool to measure the electronic structure on the surfaces of conductors. Atomic resolution and even manipulation of 
single atoms on a surface has been achieved in many laboratories [Q. Most often STMs are used in a two terminal 
configuration, the probe being one contact and the STM tip being the other. It has been shown that at zero 
temperature the two-terminal conductance between tip and surface is proportional to the local density of states 
(LDOS) on the surface at the coupling point x of the tip, v(x), and given by the Bardeen formula |2j, Gu v>sur f — 
(e 2 /h)^ir 2 vtip\t\ 2 v(x). Here, v t i V is the density of states in the tip and t is the coupling strength of the tip to the 
surface. 




FIG. 1. a) Experimental setup of an STM used to measure on a mesoscopic multiprobe conductor. The wire is biased with 
potentials V\ and V2, the applied potential at the tip is Vu p . This setup can be used to measure the average current or the 
fluctuation spectrum of the current flowing into the tip. b) Spatial variation of the current fluctuations and L(P)DOS of a 
ballistic wire with a 8 barrier at x = with transmission probability T — 0.7. The injectivity of the left contact, v(x, 1), 
(dashed line), and the LDOS f(x), (dotted line) are in units of 1/hv. The solid line gives the current fluctuation spectrum 
<(A/ tlp ) 2 ). 



It is the purpose of this work to suggest and theoretically investigate experiments which use an STM to measure 
on multiprobe conductors as depicted in Fig. la). By applying different voltages at the contacts of the conductor 
one can drive a current through the sample already without the presence of a tunneling contact. Transport through 
multiprobe samples is characterized by a conductance matrix G a p = (e 2 / 'h)T a p describing the conductances between 
each two contacts a and (3 of the sample. In the presence of a tunneling contact the overall conductance matrix of 
the sample plus tip now includes the conductances Gu P , a from a sample contact to the tip and the conductances 
G a ,tip from the tip to a sample contact. In contrast to the typical STM measurement which is characterized by one 
conductance only, we need for the arrangement of Fig. la) to determine two or more conductances. One aim of 
the present work is to investigate the densities which determine the transmission probabilities T t i P . a and T t i P . a for 
transmission from and to the tip. Once the conductance matrix is known we can use Biittiker's formula Q| which 
expresses the current flowing into the contacts a of an arbitrary phase-coherent multiprobe conductor with the help 
of its scattering matrix and the Fermi distribution functions f a (E) of the electrons in the reservoirs || , 

< 7 «> = ill [ dET aP {E){f a {E) - f f3 (E)} . (1) 

Contacts weakly coupled to a multiprobe conductor have been investigated by Engquist and Anderson using 
scattering matrices to describe the splitting of currents between conductor and tip. This discussion allowed a derivation 
of Landauer's resistance formula || without appealing to screening arguments (local charge neutrality). Only currents 



1 



but not their quantum mechanical amplitudes where considered in Ref. [Q. Subsequently, Imry Q applied the Fermi 
Golden Rule to investigate this problem. A fully phase coherent discussion of a weak coupling contact (also using 
scattering matrices) was provided by one of us j?]. Here we return to this problem. Instead of using scattering 
matrices which describe the coupling of the tip to the sample and a separate scattering matrix for the sample, we 
start from the overall scattering matrix which includes all components of the system under investigation, the sample 
and the tip, and use a Green's function technique to arrive at the results (s|. Such an approach permits readily to 
express results either in terms of scattering matrices or wave functions. 

Average current — We find that the transmission probability from one contact a into the tip is no longer 
proportional to the total LDOS but only to a part of it, a local partial density of states (LPDOS) u(x, a) p,pO|, called 
injectivity of the contact a at the coupling point x. The transmission from the tip into a contact a is proportional to 
another LPDOS v{a,x), called emissivity [p|,|i"l| of the point x into contact a. We find || 

Tup.a = 4n 2 i>ti P \t\ 2 v(x,a) and T a .u p = Air 2 v{a, x)\t\ 2 v t i P ■ (2) 

These results can be viewed as generalized Bardeen formulae for tunneling into multiprobe samples. The injectivity 
f(x, a) can be expressed with the help of the scattering states ip am which represent an incoming electron from one of 
the open channels m of contact a scattered in all channels of all contacts, v(x, a) = ^ m£a (l / 1 hv a m)\ipam{x)\ 2 . Here, 
l/hv am is the density of scattering states in channel m of contact a and v am is their Fermi velocity. The injectivity 
of a contact a thus gives the density of states at the Fermi energy and at the position x for charge carriers which 
entered the sample through contact a. Injectivity and emissivity obey the sum rule ^ Q v{x, a) — ^ Q v(ot, x) = v[x) 
so that Eqs. (^) reduce to the well known Bardeen result in the case where the conductor is only connected to one 
single electron reservoir. Whereas the LDOS v{x) is an even function of the magnetic field £>, the injectivity and 
emissivity are related to each other by the symmetry relation vb{o.,x) = V-b{x,cx). The emissivity can therefore 
be expressed with the help of the scattering states of the Hamilton operator in which the magnetic field has been 
reversed. In addition, this symmetry relation manifestly shows that the transmission probabilities, Eqs. (||), satisfy 
the Onsager-Casimir reciprocity relation ||, T a ^i P {B) — T t i P , a {~ B) . 

Using (||) in Eq. (Q) we can express the current flowing into the tip in terms of the applied potentials and the 
injectivities of the sample. At finite temperature we have 

(Iu P ) = lj dET ts (x){f tlp (E) - f eff (x)} (3) 

with the two-probe tip-to-sample transmission T ts {x) — ^T: 2 vtip\t\ 2 v{x) and the effective local distribution function 

/e//(*)=£^/^)- (4) 

This expression gives the local non-equilibrium distribution of charge carriers at the point x inside the conductor. Its 
energy dependence comes from the Fermi distribution functions and a possible energy dependence of the L(P)DOS. 
Pothier et al. Jl2[ measured an over a spatially wide range averaged non-equilibrium distribution function using a 
large superconducting tunneling contact on a diffusive wire. Eq. ([}]) has the form of the current in a two probe 
system, one probe being the tip, where the electron distribution is described by the Fermi function fu p {E) and an 
other probe where the electron distribution is given by the effective distribution function f e ff(x). This effective 
distribution function does not account for any energy relaxation of the charge carriers inside the conductor. We 
assume that electron-electron and electron-phonon interactions can be neglected for the system under consideration 
and therefore the phase of the charge carriers is conserved. However, the distribution function does contain via the 
L(P)DOS the quantum mechanical phase coherence of the electron wavefunction throughout the system. Our effective 
distribution can be used to describe the quasi-particle distribution in phase-coherent diffusive conductors, if energy 
relaxation and dephasing can be neglected. To describe transport and noise in diffusive conductors one can also use 
the semiclassical Boltzmann-equation approach, which introduces a distribution function which does not contain the 
quantum-mechanical phase-coherence but where energy relaxation processes can be modeled quite easily. However, 
the distribution function of this semiclassical approach can not be used for conductors where phase-coherence is 
essential. 

At zero temperature we can replace the Fermi functions in Eq. (|^) by step functions and get in linear response to 
the applied potentials 

(I tip )=G(x){V tip -V eff (x)} with V eff (x) = J2^^V a , (5) 
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and the conductance G(x) — (e 2 /h)T ts (x) has to be taken at the Fermi energy. A particularly interesting setup is, 
when the tip is used as a voltage probe, i. e. we demand that there is no net current flowing into the tip, (Jup) = 0. 
Similar experiments, also called scanning tunneling potentiometry, have been reported in p3[ . From Eq. (pf) we find 
that at zero temperature the voltage one has to apply at the tip to achieve the zero-current condition is exactly the 
effective voltage V e ff(x) defined in Eq. (||). The injectivities v{x,a) are determined by the equilibrium electrostatic 
potential U(x) in the sample and, therefore, also the measured effective potential V e ff(x) depends on the electrostatic 
potential. However, there is no direct relation between the measured potential and the actual electrostatic potential 
in the sample j|. 

In the following we want to evaluate the effective distribution function or voltage for a metallic diffusive wire 
and for a wire with a barrier. First, we consider a diffusive wire of length L to which two electron reservoirs are 
attached at the left end (x=0) and at the right end (x=L). The injectivity of the left contact (contact 1) and the 
right contact (contact 2) is then averaged over many different disorder configurations v(x, 1) = vq{L — x)/L and 
v(x,2) = vqx/L with the constant LDOS vq. This gives for the effective distribution function the known classical 
result f e ff(x) = f\(E){L — x)/L + fe{E)x/L. It is interesting to contrast this result, with a a fully phase-coherent, 
sample specific result. We investigate the potential measured on a ballistic one-channel wire with a barrier at x = 
leading to the transmission probability T and reflection probability R = 1 — T. Such a barrier generates strong, 
Friedel-like oscillations [0J|] in the neighborhood of the barrier. To the left of the barrier, the injectivity of the left 
contact, contact 1, and right contact, contact 2, are 

v(x,l) = + R + 2\ r Rcos(2k F x + <l))) and v(x, 2) = —T, (6) 
hv hv 

were <fi is the phase acquired by reflected particles. At zero temperature this gives the oscillating effective voltage 
V e f f (x) =V\ — T{V\ — V2) /{2 + 2\/Rcos(2kFX + (/>)}. The voltages measured to the left and to the right of a scattering 
region can be used to find the four-probe resistance of the scattering region, which is defined as the measured potential 
drop divided by the current flowing through the scattering region || . The measured resistance will in general strongly 
depend on the exact positions where the voltage is measured as the strong oscillations in V e f / (x) show. 

Current fluctuations — Until now, we were only interested in the average currents given by Eq. (|l|). Now, 
we want to go one step further and investigate the fluctuations of the current away from its average value. The 
current fluctuation spectrum can give more information about the transport properties of a conductor than can be 



drawn from pure conductance measurements. Van den Brom and van Ruitcnbcck 1 14| used combined conductance and 
current fluctuation measurements to characterize the detailed transport mechanism through few-atom gold contacts. 
Shot noise measurements have also been used to identify the fractional charge of the quasiparticles in the fractional 
quantum hall regime p"5[ . The low-frequency current fluctuation, respectively, correlation spectra are defined as the 
Fourier transform of the current-current correlator and given in terms of the scattering matrix s Q ^ of the system as 

(M a M p ) = J dt{AI a (t + to)M p (to)) = \ Y,j dETr[A Sl (a)A lS (f3)]f s (l - / 7 ) . (7) 

<5 7 

with AI a (t) = I a (t) — (I a (t)) and the current matrix Af} 1 {a) — 8 a p5 ai — stgS a7 . If the two indices a and j3 are 
equal, then this equation gives the current-fluctuation spectrum at the contact a. Otherwise, the equation gives 
the correlation spectrum of the currents at contacts a and (3. Current fluctuations at the contacts of a mesoscopic 
conductor have been measured by many groups |17|. Birk et al. f]l| used an STM to measure the shot-noise on a 



tunneling contact. Correlations of currents at two different contacts of a sample have been measured only recently by 
Henny et al. @ and Oliver et al. @. 

We calculate the current-fluctuation spectrum at the tunneling contact for the setup shown in Fig. la). As in our 
discussion of the sample to tip conductances, Eq. (^|), we start with the scattering matrix of the entire system (wire 
and tip) and develop the scattering matrix elements describing scattering from and to the tip to the lowest order in 
the coupling energy \t\ of the tip. We assume that the potential at the tip is adjusted such that the average current 
at the tip vanishes. Then, using G(x) = (e 2 / h)Air 2 vtip\t\ 2 v(x) , we find for the fluctuations of the current at the tip 

((AI Up ) 2 ) = 4 J dEG(x){l - f tip (E)}feff(x) • (8) 

We see, that the fluctuations are also determined by the effective distribution function, Eq. (^). 

At elevated temperatures the fluctuations are due to thermal noise in addition with an excess noise, called shot- 
noise, which is due to the discreteness of the charge carriers. In Eq. M) the integral over energy is from the bottom of 
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the conduction band to infinity. At a temperature T and applied potential differences AV, the relevant contribution 
to the current fluctuations comes from the integration over an energy intervall of AE w max(eA V kT) around the 
Fermi energy. If the L(P)DOS are (nearly) independent of energy in this energy range, we can evaluate the integral 
over products of Fermi functions and get 

{ (AI Up f) = 2e G WA^^ (coth ^ _ J + coth __ j | (9) 

with AV = V\ — V% the applied potential at the wire. In the limit of high temperature, kT 3> eAV, this equation 
leads to the Johnson-Nyquist formula for thermal noise, ((AI t i P ) 2 ) — 4eG(x)kT. At zero-temperature, we are dealing 
with pure shot-noise. If we restrict ourselves to small potential differences so that linear response theory is justified, 
the fluctuation spectrum is completely determined by the properties of the system (its LPDOS) at the Fermi-energy 

13, 

{ (AI Up f)=4eG( X )AV^§ (l - ^) • (10) 

v[x) V nx) J 

Using the injectivities, Eq. (^), the measurement of the current fluctuation spectrum, Eq. ([Tof), is illustrated in Fig. lb) 
for the case of a measurement on a ballistic one channel conductor with a (5-barrier at x — leading to the transmission 
probability T = 0.7. 

In conclusion, we presented generalized Bardeen formulae which for the transmission from the sample contacts to 
the tunneling tip are related to a local partial density of states, the injectivities v(x,a), and for the transmission 
from the tip to the sample contacts are related to the emissivities. We have shown that measurements of the average 
current and the fluctuations of the current at a tunneling tip which probes a phase-coherent multiprobe sample are 
determined by an effective local non-equilibrium distribution function. This distribution function is given as products 
of injectivities and the equilibrium Fermi distribution functions in the electron reservoirs. We illustrated the results 
for the case of ballistic conductors with barriers and metallic diffusive wires. 
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